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c/2 ■ Abstract 

'>'■ 
k><( \ In the framework of the reduction technique for Poisson-Nijenhuis structures, we derive 

H \ a new hierarchy of integrable lattice, whose continuum limit is the AKNS hierarchy. 

In contrast with other differential-difference versions of the AKNS system, our hierarchy is 

endowed with a canonical Poisson structure and, moreover, it admits a vector generalisation. 

We also solve the associated spectral problem and explicity contract action-angle variables 

through the r-matrix approach. 



1 Introduction 

The search for discrete integrable systems got a novel impetus in recent years: the 
quantisation of integrable PDE.s [1-4], the recent findings on integrable discrete-time 
systems [5 - 7], including the "extreme" case of cellular automata g], and even some 
results in string theory and 2 — D quantum gravity (related to the so-called "discrete 
string equations" ) are perhaps the main motivations for such a growing interest in that 



field. The results reported in this paper can be ascribed to the above line of research, 
although they belong in a sense to a more traditional approach, aiming at deriving 
integrable systems with discrete space and continuous time. 

In fact, we present and discuss an integrable differential-difference version of so- 



called AKNS hierarchy M, already mentioned in 10]. In comparison with other dis- 



cretisations |11| a, it has certain advantages and some drawbacks. The latters are 



essentially given by the non-existence of one-field reductions, unlike the model de- 



rived by Ablowitz and Ladik |ri|: such reductions are in fact admissible only in the 
continuum limit. However, this seems the price to be paid (i) to preserve the first 
hamiltonian structure of the AKNS hierarchy, namely the canonical one, and (ii) to 
allow for a vector generalisation, both such features being exhibited by our model. We 
have to mention that one equation belonging to our hierarchy can be interpreted as a 
Backlund tranformations for the continuous AKNS system, and, as such, it has been 
recently derived by Yamilov and Svinolupov |1^ . 



The paper is organised as follows. 

In Section 0, the hierarchy under scrutiny is derived by using a nowadays standard 
geometrical reduction technique for Poisson-Nijenhuis structures. 

In Section |^, the underlying direct and inverse problem is solved: for the sake of 
simplicity, only the 2x2 matrix case is considered. 

In Section ^, the r-matrix structure of the system is revealed, and, through the 
r-matrix, the action-angle variables are explicitly given in terms of the spectral data. 

In section |^, it is shown that the whole hierarchy goes into the AKNS one in a 
suitable simple continuum limit. 

In Section ^ we make a few concluding remarks and mention some interesting open 
problems. 

2 Poisson-Nijenhuis structure 



In a recent paper, P.M.Santini and one of the authors []T3|, have considered the following 
abstract linear problem: 

Ei; = {Q + \A)ij (1) 

In formula (|l|), A G C is a spectral parameter, Q, A, ip take values in an associative 
algebra A with unit element /, endowed with a trace-form < -, ■ >, i? is an algebra 



automorphism; A is assumed to be a constant (= fixed point of E) element in tlie 
algebra. 

It lias been shown in JlSl that to the linear problem (|l|) one can naturally associate 
the following two compatible Poisson tensors: 

di = LqRa — RqLa (2) 

02 = MOi (3) 

where the symbol L^ (resp.i?^;) denote left (resp. right) multiplication by x, and M is 



the hereditary recursion operator, or Nijenhuis tensor |JTJ], defined as: 

M = {LqE - Rq){LaE - RAr' (4) 

Furthermore, one has proved there that the family of vector fields: 

Kb = {Rb-Lb)Q=[Q,B]; [A,B] = (5) 

is an invariant family for Af, namely the Lie-derivative of TV along Kb vanishes: 

Cks^ = (6) 

We can then assert that to the linear problem (|I]) it is naturally associated the following 
hierarchy of (commuting) bi-hamiltonian systems: 

wr^'"^ (7) 

Out of the abstract hierarchy (J^), one can construct bi-hamiltonian lattices through 
a convenient realization of the algebra A, chosen to be the algebra of matrix valued 
sequences, approaching an arbitrary constant value as \n\ -^ oo, and of the automor- 
phism E, taken to be the "shift operator" on sequences: 

EXn = Xn+1 neZ (8) 

In particular, throughout this paper, we shall consider (A^ + 1) x (A^ + 1) real matrices. 
Accordingly, the "field variable" Q will be parametrized as follows: 



where p is a scalar, < q\ is the 1 x A^ matrix (gi, ^q^) ("row vector"), \r > is the 

A^ X 1 matrix (ri, ,r/vr) ("column vector") and Q is an A^ x A^ matrix. 

The constant matrix A will be chosen as: 



^= ,., .' (10) 




where < 0| (|0 >) is the null row (column) vector, and is the null A^ x A^matrix. 

The sequence {Qn}n£Z of matrices of type (|^) is assumed to fulfil the boundary 
condition: 

lim Qn=f ° 7') (11) 

where / denotes the N x N identity matrix. 

Let X C ^ be the subalgebra of A, consisting of matrix- valued sequences obeying 
homogeneous boundary conditions. Once equipped with the Lie-product given by the 
point-wise commutator: 

I becomes a Lie-algebra and moreover an ideal of A, and can be identified with its 
dual through the bilinear form: 

+ 00 +00 

{X,Y)= Yl <Xr,,Yr,>=Tr Y. X^Yn (12) 

n=—oo n=—oo 

Our configuration space M will be the affine hyperplane to I of matrix- valued sequences 
obeying (pJ]); accordingly, its generic point will be again denoted by Q. 

Its tangent bundle T{M) will then be the set {(Q, Kq) : Q e M,Kq El} while 
the points of its cotangent bundle T*{M) will be parametrized as {(Q,7q) : Q G 
M,7qGJ*~X}. 

In the following, for the sake of simplicity, elements of T{M) (resp. T*{M)) will be 
shortly denoted by Kq (resp. -yg). In the present concrete realisation of the abstract 
setting introduced in |jl3|, the Poisson tensors 9i, 62 defined by (^,H), as linear maps 
form T{M) and T*{M), are in fact linear operators on I, and the same is true for the 
Nijenhuis tensor A/", given by (^, which is a linear map form T{M) into itself. 

The main result of present Section is contained in the following: 

Theorem 1 The Nijenhuis tensor M reduces by restriction on the submanifold Im 61 . 



As a consequence, the abstract bi-hamiltonian hierarchy (0) restricts to a bi-hamiltonian 
hierarchy of evolution equations on a one- dimensional lattice. 

The theorem (|l|) will be proved in a constructive way, which will then also yield 
the concrete form of the restricted Poisson and Nijenhuis tensors. 

Let us start by evaluating the image of 9i: 

Imd, = {KQ = eaQ;iQeT} (13) 

Let us choose for Kg and jq the parametrisation induced by (^): 

rr I Kp <Kg\\ ( -fp <lr\\ ...^ 

One immediately gets, for Kq G Im6i, i.e. for K of the form: 

Kq = QiqA - A^qQ (for some 7q) (15) 

the necessary condition: 

Kq = (16) 

Hence, Q must be a constant N x N matrix. In view of the boundary condition (p!TD, 
we have: 

Q = i (17) 

In turn, formula (p!7| ) implies: 

Kp =< Kq\r > + < q\K,. > 

which entails: 

p = qr + c (c = const.) (18) 

Summarizing, we can assert that: 



1. under conditions (|10D, ([llD , Im6i is the submanifold of X consisting of vector 
fields Kq such that: 

Kq = 6 ; Kp=< Kq\r > + < q\Kr > (19) 

2. the set 5* C M, such that T{S) = ImOi , given by: 

S = {Q:p=<q\r>,Q = i} (20) 

is a characteristic leaf of 9i. 



We are now ready to show that TV reduces by restriction on S, namely that T{S) = ImOi is an 
invariant submanifold for TV. To this aim, let us introduce the auxiliary vector field ^q through the 
formula: 

Kq = {EiQ)A - ACq (21) 



whence it follows, on T{M): 
Parametrizing S,q as 

we get from (pl|) 



K'Q^MKQ = {Ei^)Q-QS,Q (22) 



i. - l^ <f I (.3, 



6 = {E-l)-^Kj, (24) 

<e,| = -<K,\ (25) 

\ir > = E-'^\Kr > (26) 



while the matrix S, stays undetermined. 

On the other hand, eq.(|2^) yields, for Q ^ S: 



K; = < q\r > {E - l)Kp+ < E^,\r > - < q\^r > (27) 

<m - < q\{E^a) + (E- < q\r >) < ^g\- < q\i (28) 

\K> = i< q\r > E ~ mr > -Ur > HEi)\r > (29) 

By imposing the reduction condition: K'q — 0, one gets: 

i=iE-l)-\-\E^r><q\ + \r><C,\) 

or, using (||,|^: 

e - -{E - ir\+\Kr X q\ + \r >< K,]) (30) 

Then, to establish that S is an invariant submanifold for A/", we have just to show that: 

Kp =< Kg\r > + < q\Kr > (31) 

entails: 

K; ^<K'^\r> + < q\K', > (32) 

and this can be seen by a direct straightforward (although tedious) calculation. □ 



Hence, the theorem ([1|) is proved and moreover we can give the esphcit form of 
M\s, more precisely we have: 



\K> 



i<q\r>-E)\Kg>+(^^^{<Kq\r>+<q\Kr>)+^^i\r><Kq\ + \Kr><q\))\g> 
{<g\r>-E-'')\Kr>-{^{<Kg\r>+<q\Kr>)+^(\r><Kg\ + \Kr><q\))\r> 



(33) 



The exphcit form of ^il^ can be easily obtained by noting that the points of S 
(resp. T{S)) are completely determined once the 2N vector (< q\,< r|)* ( resp. 
(< Kq\,< KrlY) are given. 

Hence points of T*{S) will be fully characterized by the value of the 2N vector 
{< Pq\,< Pr\) defined by the duality condition: 



Tr J2 lQ,nKQ,n = E <^q^Kq„> + < /3rjKr„ > 



(34) 



which entails: 



< f3q\ =< 7^1 - 7p < r| < f3r\ =< 7g| - % < q\ 

By direct calculation, one can check that the relation: 



once restricted to S, implies: 



\Kq> 
\Kr> 



-/ 
/ 



1/5. > 
1/3, > 



= J 



\Pq> 
\/3r> 



(35) 



(36) 



whence it follows that 6*115 is the canonical Poisson tensor or, in other words, that the 
fields variables \q >, \r > are endowed with canonical Poisson brackets: 



{(lj{n),qk{m)}i = {rj{n),rk{m)}i 







(37) 



It is worthwile to notice that the restricted Poisson tensor 6*215 = -f^ls^ils is on the 
other hand nonlocal and has a rather cumbersome form: its non-locality is intimately 
related with the non degeneracy of 6*115. 

According with the general theory of bi-hamiltonian system, we have still to show 
that the restricted vector fields: 



kB = [Q\s,B] {[B,A]=0) 



are invariant vector fields for N'ls'- but tliis follows by the invariant nature of the Lie 
derivative, provided that Kb belongs to T{S). 

Actually, the commutativity condition [B, A] = implies for matrix B the form: 



B 



b <0| 
10 > B 



so that, for any Q G 5, we have: 

Kb = 



<q\{b-B) 

{B-b)\r> 



(3^ 



So, the family of starting commuting symmetries can be written in terms of the 

2A^-vectors: 

' C'\q> 

-C\r> 






(39) 



C = b — B being an arbitrary N x N matrix. 

However, we should notice that the resulting evolution equations: 



\q> 

\r > 



ms 






(40) 



will be in general non-local. Local equations will be get whenever C = cI. 

Let us now give some concrete examples of local evolution equations associated 
with the linear problem ([l|), corresponding to the choice = 1: 



1. {I = 0) 



or, in components: 



\q> 
\r > 



to 



\q> 
-\r > 



n\s 



\r > 

\q> 



dq^^ 
dto 

with hamiltonian density: 



?i 



^ = -r,- (j = l,....N) 



N 



ho 



< q\r >= — ^ qji 



j=i 



(41) 



(42) 



(43) 



2. (/ = 1) 

|g > \ / (< Q.V' > -E)\q > 



1 ._,._.^,._. , (44) 



\r > I \ (E — < q\r >)\r > 



with hamiltonian density: 

1 



hi = -^ < q\r >^ + < Eq\r > (45) 



3. (/ = 2) 



q> \ I -E'^\g>+E{<q\r>\q>)+{<Eq\r>+<q\E-'^r>-<q\r>^)\q> 

(46) 



|r> / \ \E-^-<Eq\r>+{<q\r>)^-<q\E'^r>j\r>-[<q\r>+<E-^q\E~'^r>j\E-'^r> 



with hamiltonian density: 

/i2 = — < q\r >^ + < q\r > (< g|E"V > + < Eq\r >)- < Eq\E'^r > (47) 
o 

4. (/ = -l) 

k>] ^ J \q>\f-E-^\q>/il^<E-^q\r»\ 

\r>)^ \-\r> ) \ E\r> /il+<q\Er>) J ^ ' 

with hamiltonian density: 

/i_i = ln(l+ < q\Er >) (49) 



We will now show that the above geometric reduction procedure is perfectly equivalent 
to the peharps more familiar technique, based on the so-called discrete zero-curvature 
condition: 

Ut + UV - {EV)U = (50) 

obtained by enforcing compatibility of linear problem (|l|) with the auxiliary linear 
problem: 

^ = V{X)ij (51) 

In our case: 

U = Q + XA (52) 



where A is given by (0) and Q belongs to manifold S (^0]), i.e.: 

Q=\ <^l'> <.^l 
\r > I 



To extract from (g) the hierarchy (|40|) we start by considering the "stationary equation" 
associated with (pDI), namely: 



UW - {EW)U = (53) 

Parametrizing W as: 

^=f " ^^"'1 (54) 

one gets 

# = (E - l)-i[|r >< u\ - \Ev >< q\] + W{X) (55) 

w = -trW = {E- 1)-% >< m| - \Ev >< q\] - tr W{\) (56) 

W{X) being an arbitrary constant (i.e., field independent) matrix, and the following 
"eigenvalue equations" for < u|, |t> >: 

X<u\ = {E- < q\r >) < u\+ < q\{Ew - W) (57) 

XE\v >= {E'^- < q\r >)\Ev > +{w - EW)\r > (58) 



with w, W given by (|55|) , (pq) ■ 

Eqs. (|57| , p8D can be solved recursively assuming W, and thus < u\, \v >, to be 
given by the following Laurent series in A: 

oo 

J^iy^'^U-'^ (59) 

k=0 

and requiring W{\) in (|55|) , (|56|) , to be in fact A-independent. 
One gets: 

< M^'^+i)] = {E- < q\r >) < m(^)|+ < q\{Ew^''^ - W^''^) (60) 

^|^(fc+i) >= (^-i_ < g|^ >)|Et;W > +(u;('^) - E#('=))|r > (61) 

On the other hand, if PF is a solutions of (|53D , the same is true of course for X'^W, 
where k is any positive integer. 
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Given a two-sided Laurent series on the unit circle: 

+00 
/(A) = E fs>^' > 

s=— 00 

we shall denote, as usual, by /+(A) (resp. /-(A)) the part containing only non negative 
(resp. strictly negative) powers of A. Therefore, eq. (^) can be rewritten as: 

U{X''W)+ - {X''EW)+U = -U{X''W)- + {X''EW)-U (62) 

But now, as U is linear in A, the l.h.s. of (^) cannot contain any negative power of A, 



while its r.h.s. cannot contain any strictly positive one. Hence both sides of (|62|) are 
A-independent (order "zero" in A), and we have: 

l.h.s. = QW^''^ - {EW^''^)Q (63) 

r.h.s. = -AW'^''^^^ - {EW^''+^^)A (64) 



As far as its A-dependence is concerned, (|63|) (and of course ^4]) is then compatible 
with Ut. On the other hand it may be readily checked that it belongs indeed to the 
manifold T{S), defined in eqs. (p!9|-pO|). Hence, we can assert that the hierarchy of 
evolution equations associated to (|l]), (^Tj) is given by: 

f)TT 

^ = -QiyC^) + (Eiy('=))Q = Aiy^'^+i) - (EW^''+^^)A (65) 

Otk 

They correspond to the following choice for the matrix V appearing in formulas (^Tp, 

my- 

V = %) = iX^W)+ (66) 

and clearly coincide with (^), by choosing — C = (trW)I + W{see rqs. (|39D, (|55| ) ,(|56| ) ). 
In particular, C = I '^ W = —j^ril. Indeed, looking at (^), (p3D , we see that W is 
the generating function of the fields ^q. 

3 Direct and inverse problem 

In this Section, we outline the solution of the direct and inverse problem associated to 
(0); for simplicity, we restrict considerations to the 2x2 matrix case, when we have: 

U„iX) = Q„. + XA=\ ^ + ^""" ^; I (67) 
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As explained in section 2, {qn}n£Z, {rn}n£Z are real valued sequences vanishing at ±00 

lim Qn = lim r„ = (68) 

Inj— >oo |n|— >oo 

For the linear problem: 

V'n+l = Untpn (69) 

we can naturally define the transfer matrix: 

Tn,mW = Un-1 Um {u > m) (70) 

-1 

2+1 

such that: 



T„,^(A) = U-'U-l, f/,;ii [n < m) (71) 



We can then introduce the ^^Jost matrices" : 

T-(A)= lim T„,„(A)E^(A) = (0„(A),0„(A)) (73) 

m— >— oo 

T+(A)= lim T„,^(A)E„(A) = (0„(A),7A„(A)) (74) 

m — ^+oo 



where: 



En{X) = I J (75) 



is a fundamental matrix solution of the asymptotic (or "undressed" ) problem: 

En^lW = ( 1 ) ^«(^) ^^6) 

and ifn, <pn, ^n; '^n are 2-column vector solutions of (|I|), the ^^Jost solutions" . 

Cleary, the asymptotic solution -En (A) ([75|) is bounded on the unit circle |A| = 1, 
which will then be the continuous spectrum of (|l]), (|68D . 

On the unit circle, the monodromy matrix is then defined as: 

T(A)=T+(A)[T„(A)]-^ (77) 

In the following, we shall call ^^ spectral parameters" the elements of monodromy matrix: 
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As det Un = A, we have: 

det T„,^(A) = A"-'", det T±(A) = A" (79) 

so that: 

det r(A) = 1 (80) 

Formulas (^) imphes that </?„, (pn and ipn, (pn are two pairs of independent vector 
solutions of (0) on the unit circle (|A| = 1), while formulas ([78|) , (|80|) entail, on the 
unit circle: 

a(A)a(A) = 1 + b{X)b{X) (81) 



Direct problem 

The direct problem amounts to determine the monodromy matrix ( [78| ) once the 
fields {gn}, {'"n} are given. 

To this aim, it is convenient to introduce the normalized vector sequences: 

X = A->„; ^ = A->„ (82) 

such that: 

lim Xn(A) = ( M = lim ¥^„(A) (83) 

Somewhat loosely, in the following Xn, V'n will be denoted as "Jost solution" as well. 

The spectral parameters are naturally defined as Wronskians of independent vector 
solutions; namely we have: 

a{\) = Wi^,i;) 

~a{X) = W{4>,x) .g^. 

b{\) = X^Wix,v) ^ ' 

6(A) = A-"iy(0,V) 

where W{a, b) is the determinant of the matrix whose columns are the 2- vectors a and 
b. 

It easily seen that the vector sequences Xn, i^n, ^n-, <^n satisfy the following 
"discrete integral" equations: 
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V'n(A) 
0n(A) 



i)\?J 1 



1 





1 



k=n 
n-l I ^-1 



E 







A 



k—n 



fc=— oo 
n— 1 / \n~k—l n 

1 



E 



rfc 

QkTk qk 

Tk 

QkTk Qk 



^fc(A) 
V5fc(A) 
0fc(A) 



The analyticity properties of the Jost solutions are summarized by the following: 
Theorem 2 If the sequences {g„}, {r„} are such that: 

lim n^g„ = lim n^r„ = 

|?i|— >oo |n|— >oo 



(86) 
(87) 
(88) 



^9) 



then ipn, fn OLfG analytic functions of X in the domain |A| > 1 and are continuously 
differentiable for \X\ > 1; analogously, Xn, 'fn are analytic functions of X for \X\ < 1 
and continuously differentiable for \X\ < 1. 

We outline the proof of theorem (|^) for (y9„. 
First of all, we equip C'^ with the norm: 

||a;|| = max \xk\ 

k=l...N 

so that linear transformations in C'^ are naturally equipped with the norm: 



N 



1X11 = 



max > 

3 = 1. ..N ^ 



J = l 



Then, the following inequality holds for (87): 

A-i \ 



A 



k—n 



Ui, 



< 



\\Ui,k\\ (|A| > 1) 



where 



Uuk = 



qkTk qk 

Tk 



So, writing the Neumann series solutions of ( p5|) 



M.,^ £<-«"(;) 



where: 



f(^) 



1 
1 



(90) 



(91) 



(92) 
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^i'^A)- Y, G{nM) E G(fci,fe) Y. G{h^i,ki) (l>{)) (93) 

fci>n fe2>fe2 k2>k2 



A-1 
A* 



GKfc)-( . ,.-„ ]Ui,k (94) 



The inequality (^0|) implies, by iteration: 

\\FiHm<^(Y.\\Ui,k\\] l>0 (95) 

Therefore: 



\k>7i 



||(p(A)||^sup||</p„(A)||<exp7 (96) 

nez 

where: 

nGZ 

Hence, if 7 < 00, <y9n(A) is analytic for |A| > 1. On the other hand: 

||C/i,„|| - niax{|(7„|[l + |r„|] , |r„|} < |9„| + |r„| + 1/2[|9„|2 + |^^|2] 
so that the existence of 7 is guaranteed whenever the sequences {(?«}, {''«} belong to h. 

d\ 



A similar procedure leads to the following result for ^^ ■ 



dX 
which holds, with suitable coefficients a and 6, whenever 



<a+|n|6 |A| > 1 (98) 



J2i^ + \n\)\\U^,,\\<oo (99) 



nez 

,-2 



Hence, provided qn and r„ vanish faster than n as \n\ — > cx), (fniX) is continuously differentiable 
with respect to A for |A| > 1. □ 

We can thus assert that, whenever {g„}, {r„} belong to h, the diagonal entries 
of the monodromy matrix a(A) and a{X) (3.14) are analytic respectively outside and 
inside the unit circle. Morever, due to their asymptotic behaviour in A: 

lim a(A) = lim det{ip,ilj) = det I 1=1 (100) 

|AHoo |A|^oo \ 1 



lima(A)= lim det(^,x) = if (1 + r,g,-i)rfet ( ^""^ ^ 1= (101) 

j=-oo 



- n (i+o?.-i) (102) 



j=-oo 
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both a and a have at most a finite number of zeros, say A^ and A^ respectively, in their 

\,}JL, and {A,}^^ 



analyticity domains. These zeroes will be denoted as {\j}^^i and {\j}^^^ and will be 



assumed to be simple. 

If in addition the stronger condition (p9| ) is satisfied, then the entries of the mon- 
odromy matrix ( ffSD are Holder continuous on the unit circle (see again (|8^). Morever, 
the analyticity properties of a(A) and a(A) imply that the scalar Riemann problem on 
the unit circle (|8H) can be solved through the formulas: 



a(A) 



Y\U^ 



A,) 



n^li(A-A, 



A'^exp 



1 r ln(l + 6(;z)6( 
2'Ki J\z\=i 



alnz 



z-X 



-dz 



(103) 



for |A| > 1 



a(A) 



nf=i(A-A, ) 
nf=i(A-A, 



exp 



1 r hill + b{z)b{z)j — alnz 
2'Ki J\z\=i 



z-\ 



dz 



(104) 



for lAI < 1 



where a is the index of the Riemann problem, i.e. the variation of Arg{\ + 6(A)6(A)) 
after a cycle; condition ( |100|) clearly implies a = N — N. 

In the following, we shall assume the index a to be zero so that N = N: this is 
certainly true in the refiectionless case (6 = 6 = 0). 



Inverse problem 

The inverse problem amounts to reconstruct the sequences {g„}, {^n}, once the 
monodromy matrix is given. 

To solve it, we rewrite eq.([77|) in terms of the appropriate column vector solutions: 



^n , , ^\n~ 

— = V'n + -A"Xn 

a a 



V^n ~ 6 

— = Xn + -A "V'. 

a a 



(105) 
(106) 

(107) 



Under the above analiticity conditions on 0, ip, ip, x, eqs. ( |105| , |107|) are two vec- 
tor Riemann- Hubert problems, which can be solved (i.e. reduced to singular integral 
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equations) through the well-known Plenielj formulas. One gets: 

\l J 2m J C - A ,^^, (Afc - A) 

^"^'^ = 1, j -^r^^^T^^'^^S (A. -A) ^^°^) 

where we have set: 

p(A) = 6(A)/a(A); p(A) = 6(A)/S(A) (110) 

7k = h/a\\k); lk = hk/a'{\k) (HI) 



In formulas ( |111| ) h^ ( resp.6fc) are the ratios between cf) and ip (resp. (j) and ^) at 



A = Afc (resp. A = A^). 

Of course, in the reflectionless case formulas ( |108yi09| ) yield a system of 2N linear 
algebraic equations for '?/'n(Afc), Xn(AA;)- 

Finally, as usual, the sequences {g„}, {r„} are given in terms of the leading terms 
of the asymptotic behaviour of the vector solutions tpn-, Xn- For instance, taking into 
account that: 



^n=i *! 1+A-M ^;i) i+o(A-^ 






we have: 



Qn = -^K 



To find the time-evolution of the spectral data corresponding to eq. (0), (0), we 
notice that, if the matrix V appearing in the ausiliary linear problem (^) is given by 
(|6^), the monodromy matrix undergoes the time evolution: 

^^^^^-[%)(A),T(A)] (112) 



where: 



dt, 



%)(A) = lim %) .(A) = A'^ I "^" ^ ° I (113) 
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Restricting considerations to the 2x2 case, we can thus write: 

%)(A)=cAV3 

where c is an arbitrary scalar constant and as is the usual Pauli matrix. 
Consequently we have: 

da da 

dtk dtk 
db 



(114) 







dh 



-2c\% 



dtk 



(115) 
(116) 



of course, eqs. ( |115| ) imply that A^, A^ are constant in time for any equation of the 
hierarchy ( "isospectral deformation"), while the normalization coefficient 7^, 7r evolve 
according to equations: 



dtk 



-2(A,,)^. 



kz. 



dt, 



2(A,,)''7. 



;il7) 



4 r-matrix and action-angle variables 

We have seen in sec. ^ (formulas (p6D,(p7D) that our hierarchy of discrete evolution 
equations are hamiltonian with respect to the canonical Poisson tensor or, in other 
words, that the fields variables \q >, \r > are endowed with the canonical Poisson 
bracket. 



As consequence, we have the "ultra-local" Poisson bracket relation |T5 
{Unix) ? UUf^)} = [r(A, /x), f/„(A) ® UUf^)]6n,m 



118) 



Restricting again considerations to 2 x 2 matrices f/„, we have for 4x4 r-matrix r(A, /i) 

the formula: 

/ 1 \ 



r(A,/i) =r(A-/i) 



1 



A — /i 



10 
10 

V 1 ; 



;ii9) 



i.e. the same expression as for the NLS case. 

From (|118| ) one easily gets the analogous relation for the transfer matrix, and finally 
the Poisson-bracket relation for the monodromy matrix, that reads: 



{T(A) ? T(/x)} = r+(A, /i) [T(A) ® T{f,)] - [T(A) ® T(/i)]r-(A, /.) 



(120) 
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where: 



that is: 



r^{X, /i) = hm {E''^{n, A) ® E'^{n, /i)r(A, fi)E{n, A) ® E{n, /x) 



/ -±- 



\ 

TiTr6{X-fi) 

±i7rS{\-fx) 

V ^y 



(121) 



In formula (|121|) , j^ denotes the principal- value distribution, and we have used the 

WA)" 



X—fj. 

distribution formula: 



lim 



^in5{X — [I) 



)i-»±oo A — /i 

Eq. (|120|) implies the following Poisson-brackets for spectral parameters: 



{a(A),a(/i)} 

WA),%)} 
{6(A), %)} 

{a(A),%)} 



{a(A),a(/i)} = {a(A),a(/i)} 
{6(A),%)} = 







(122) 



(123) 
(124) 
(125) 



-{%),«(A)} 



X — ^ 
We can thus construct canonical spectral variables, given by: 

a(A) = lna(A)a(A) /5(A) = ^ In |^ 



a(A)6(/i) — i7r5(A — /i) b{fi)a{fi) (126) 



A, 
A, 



obeying the Poisson bracket realtions: 



Ani 

Vj =ln7j 



(127) 



{a(A),/3(/i)} = 5(A-/i) 



(128) 



(all other P.B. being identically zero). 

As usual, in terms of the spectral variables, the continuous and the discrete spec- 
trum contributions are separated. Moreover, a(A), {Aj}, {Aj} are constant along any 
flow of the hierarchy, while /3(A), {vj\^ {vj] evolve linearly. 

The "action" variable a(A) is deflned only on the unit circle. However, it can 
be expressed as a sum of two functions having a single-valued analytic branch for 
|A| > maxj |Aj| and |A| < miuj |Aj| respectively, uniquely deflned by the asymptotics 
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( |100| ),( p!0^ ). Such branches will be denoted as Ina(A) and Ina(A) respectively. Then 
formulas (|103|) , ([T0^ imply for such branches the following power series expansions: 



lna(A) = ^A-'=Jfc 



fc=i 



k=0 



where: 



Jk 



N 



Jk 



2m 



27ri J\z\=i 
1 



z'' ^a\z 






N 



\z\=l Z^ 



fc+1 



a\z 



M^ + ErlA-'-A-'^) (MO) 



-,k 



N 



1/1 
Jo = — : f -a{z)dz + ^(InAj — InAj) 



(129) 
(130) 

(131) 
(132) 



(A; = 0) 



In the following, we will show that the evolution equations ( |65D are generated by the 
hamiltonians J^, i.e.: 



dU 



{Jk.Un{X)} 



Indeed, from the Poisson-bracket relation (|118D , we obtain: 

1 



{trT(/i), ?7„(A)} 



A — /i 



(K+i(/t^)(AA + g) - {\A + Q)Vn{^i)) 



(133) 



(134) 



1 



Vn+l{^Ji)A-AVn{^l) + - — iyn+MU{n,ii)-u{n,^i)vM) 



X — fi 
where K, is defined in terms of the Jost solutions: 



Vnifi) 



^2,n(/U)v5l,n(/U) -V^l,„(/i)v5l,n(/i) 



^2,n(/i)v^2,n(/i) -^l,n(yU)v'2,n(/i) 

It can be immediately seen that 



-X2,n(At)01,n(/U) Xl,„(yU)0i_„(/i) 

-X2,n(/i)02,n(yU) Xl,n(yU)02,n(/i) 

(135) 



(K+i(/i)f/(n,/i) - U{n,fi)Vn{fi)) = 



Then formula (134) entails: 



{a(/i),[/„(A)} = K!:JA-AVW 
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where V"^ are of course the projections of V outside and inside the unit circle. Hence: 

{lna(/x), Un} = WltM - ^Wi+^ (136) 

{lna(/x), [/„} = l?i;U - AWi-^ (137) 



where 

,-1 



a ^V(+) 



«"- - 1 :_.;,_, (138) 

are again analytic functions of A for |A| > maxj \\j\ and |A| < minj |Aj| respectively; 
they obey the asymptotic conditions: 

lim lyW = I ^ ° 1 (139) 

|A|^oo \ 

lini ly(-) = 1 ( ^"^"-1 ^"-1 ] (140) 

|A|-o 1 + r„g„_i \^ r„ 1 y 

lim iy(+) = ( ^ ° 1 (141) 

lim ly(-) = I ° ° I (142) 

n^-oo \^ -1 y 

Formulas ([T36|),(p7D yield: 

{J,, Un} = wiXf^A - Aiy WC^) (143) 

{J,, Un} = W^-f^A - AWt^^'^ (144) 



On the other hand, W''^\ defined by (|138|) and W, defined by (p3|)(^9l), obey the same 
stationary equation ( |SBD and have the same analyticity properties with respect to A. 
Moreover, their asymptotic behaviours differ essentially just for a constant multiple of 
the identity matrix, which plays no role in the recursion relation: hence formula (|143|) 
coincides with (|65|). 
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We end this section by noting that formula ( |118| ) does hold even in the (A^ + 1) x 
{N + 1) matrix case. For this general spectral problem, the r-matrix is given by: 



/ Rl 



X — fi 



<+i \ 



\Rl 



yN+l 



N+l 



where 



{R„ 



jij OiiOfYij 



^N+1 / 



(z,j = l,....iV + l) 



(145) 



(146) 



5 Continuum limit 

Theorem 3 In the continuum limit: h —>■ 0, n —>■ 0, x = nh finite, 
with the rescalings: 

q ^ hq , r —>■ hr 

we have : 

hm = A 

h^O h 

where A is the recursion operator of the vector AKNS hierarchy. 
In fact, from eg . (|33|) : 

K'q={<q\r>-E)\Kq>+[-^{<Kg\r>+<q\Kr>)+^(\r><Kg\ + \Kr><q\))\q> 
K'r={«i\r>-E'^)\Kr>-{^{<Kq\r>+<q\Kr>)+^{\r><Kg\ + \Kr><q\))\r> 

by noting that: 

E = e^^^ = 1 + hd^ + 0{h^) 



^_i ^ ^_,,9, 



l-hd^ + 0{h^ 



and thus 



(147) 
(148) 



(149) 



^, = W + 0{1) ; ^ = 19,-1 + 0(1) 



(150) 
(151) 

(152) 
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we obtain: 



K'q=-Kq+h[-d^Kq+d^-\<Kq\r>+<q\Kr>+\r><Kg\ + \Kq><q\)]\q>+0(h'^) 



K'r=-Kr+h[d^Kr~d:,-'^{<Kg\r>+<q\Kr>+\r><Kg\ + \Kg><q\)]\r>+0{h'^) 



(153) 



so that: 



lim 



^f + I 



h^O h 



-dx+dx^^ i<-\r>+\r><-\j\q> -d^"^ (<q\->+\-><q\]\q> 
-dx''^(<-\r>+\r><-\)\r> dx-dx^'^ (<q\->+\-><q\)\r> 



(154) 



which is the recursion operator for the vector AKNS case \W\. In particular, in the 2x2 
matrix case, we get the famihar recursion operator of the standard AKNS hierarchy: 



,. Af + I 
hm 



-2rr^r- d,-2rC^q- 



h^O h 

For instance, by taking (=^^) we get the equations: 

qt = qxx - "^q^r 

< 

. rt = -Txx + 2r2g 



(155) 



(156) 



It is worth noting that the continuous equation ( |156D can be a; so obtained by taking 
the continuum limit of a suitable linear combination of the / = —1, / = 0, / = 1 flows 
(eqs. (H), O), H) ). 



6 Concluding remarks 

We would like to stress here that the system presented in this paper is, to our Knowl- 
edge, the only integrable discrete version of AKNS hierarchy that keeps the canonical 
Poisson structure of the continuum model and admits a natural vector generalisation. 

As further developments of the research reported in this paper, we mention the 
derivation of Backlund tranformations for our lattice equations, to be considered as 
integrable fully discrete two-dimensional lattices, and the formulation of a proper quan- 
tum version in the framework of the Quantum Inverse Method. 

Research has already started on both the above issues, with encouraging prelimi- 
nary results. 
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